Electronically excited cold ion crystals 
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The laser excitation of an ion crystal to high lying and long-lived electronic states is a genuine 
many-body process even if in fact only a single ion is excited. This is a direct manifestation of the 
strong coupling between internal and external dynamics and becomes most apparent in the vicinity 
of a structural phase transition. Here we show that utilizing highly excited states offers a new 
approach to the coherent manipulation of ion crystals. This permits the study of phenomena which 
rely on a strong coupling between electronic and vibrational dynamics and opens up a route towards 
the quantum simulation of molecular processes in a Paul trap. 
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Cold crystals of trapped ions are currently among the 
most precisely controllable quantum systems and are 
widely used, e.g. in quantum computation [l] [2] , opti- 
cal frequency clocks |3l H] as well as for the simulation of 
quantum magnets [SHZ], relativistic quantum mechanics 
[8 and open quantum systems [9]. All these applica- 
tions are based on the laser excitation of low lying ex- 
cited states and/or on a state-dependent coupling be- 
tween these electronic levels and vibrational modes of 
the crystal. A different regime is entered when instead 
ions are excited to high lying Rydberg states [10 . This 
promises applications ranging from fast quantum gates to 
the simulation of coherent spin dynamics and excitation 
transfer - all relying on the large dipolar interaction be- 
tween excited ions. The plethora of newly emerging pos- 
sibilities is impressively demonstrated by recent break- 
throughs in the experimental control of neutral Rydberg 
atoms demonstrating quantum gates [Til IE] and collec- 
tive excitation dynamics [13]. 

In this work we show that the laser excitation of a 
cold ion crystal to electronic Rydberg states leads to gen- 
uinely novel features which go beyond the mere extension 
of neutral Rydberg physics to excited ions. Even a sin- 
gle localized excitation exhibits a pronounced many-body 
character which is rooted in the strong coupling between 
highly excited electronic levels and vibrational modes. 
This induces state-dependent structural changes in the 
ion crystal, reminiscent of a configurational change in an 
excited molecule. Lifetimes of ionic Rydberg states on 
the order of ~ 100 /is grant long coherence times for this 
'vibronic dynamics'. This highlights novel routes for the 
coherent manipulation of cold ion crystals and paves the 
way for the study of molecular phenomena with ions in 
the well-insulated environment of a Paul trap. 

The fascinating physics emerging from couplings be- 
tween Rydberg levels and external motion has recently 
been theoretically studied for the case of neutral atoms 
by Wiister et al. [HI [15] . Here entanglement transfer as 
well as the intriguing dynamics near conical intersections 
of molecular potential surfaces were investigated. Both 




FIG. 1. Rydberg excitation of a three-ion crystal near a struc- 
tural phase transition, (a) The trap parameters are chosen 
such that ground state ions (blue) form a zigzag crystal. The 
localized Rydberg excitation of the central ion (red) effectu- 
ates a change of the structural configuration, (b) Electronic 
level scheme. For the Rydberg excitation of a "^^Ca^ ion 
from the ground state 45*1/2 the ion first is transferred to 
the metastable state 3Ds/2 via the intermediate level 4Pi/2. 
A laser of Al ~ 122 nm or a coherent three-photon excitation 
is used subsequently to populate the nPi/2 state, (c) Sketch 
of the electronic potential surface of the electronically excited 
ion crystal and a crystal with all ions in the ground state. The 
X- and y-axis show the relative displacement and the heights 
of the two potential surfaces which can be approximated as 
harmonic oscillators (here with frequencies uuq and z/q)- The 
large mutual displacement |Ko| of the two surfaces, separated 
by the energy AE, indicates that the Rydberg excitation is 
accompanied by strong mechanical forces. This has profound 
implications for the laser excitation behavior of even a single 
ion of the crystal. 

effects rely on delocalized excitations that are established 
by the dipole-dipole interaction among Rydberg states. 
While in ion crystals such a scenario can be in principle 
equally achieved [10 , the features discussed in this work 
rely on a different mechanism. Here the large polarizabil- 
ity ~ of a single ion excited to a Rydberg state with 
principle quantum number n leads to a drastic change of 
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its local trapping potential. This affects the entire crys- 
tal as changes in the position of one ion are immediately 
transmitted by the long-ranged Coulomb force. 

The phenomenon becomes most striking in a parame- 
ter regime where the ion crystal is close to a structural 
phase transition. In a linear Paul trap such transition is 
controlled by the ratio of the axial uj^ and longitudinal 
ujp trap frequencies. Depending on the so-called trapping 
anisotropy A = {oJz/ojp)'^^ the ions in their vibrational 
ground state form either a linear or a zigzag crystal . 
A transition between these two configurations takes place 
at a critical trapping anisotropy Ac ^ 2.53A^~^ ''^ [17] 
and has been experimentally observed in crystals with 
N = 3.. .19 ions [T8]. 

We are focussing here on the smallest possible system 
in which this transition is observable - a three-ion crystal 
of '^OCa+. As depicted in Fig. ^ we start in a zigzag 
configuration from where we excite the middle ion to the 
nPi/2 Rydberg state. Technically this is done by a stan- 
dard two-step excitation to the metastable 31^3/2 state 
(via the state 4P1/2) followed by a coherent excitation to 
the Rydberg level nPi/2 (see Fig. lip). The last step can 
be undertaken by either a three-pnoton excitation or by 
employing a coherent source of vacuum ultraviolet light 
with a wavelength near 122 nm as discussed in Ref. [19]. 
The excited ion will experience a state dependent - in 
general tighter - confinement 0Jp{n) which in a qualitative 
picture leads to a reduction of the trap anisotropy A{n) 
below Ac As a consequence, the electronically excited 
crystal prefers a linear arrangement of the ions. The as- 
sociated strong mechanical effects, illustrated by the far 
displaced potential surfaces depicted in Fig. will be- 
come apparent on the spectroscopic properties of a cold 
ion crystal. 

In the following we provide a more thorough quantita- 
tive study of this simplified picture. We consider a linear 
Paul trap with the electric potential <I>(R, t) = R-G-R^. 
Here, G is a 3 x 3 diagonal matrix with elements Gxx = 
acosftt — /3{l-\-e)^ Gyy = — a cos ftt — /3{1 — e)^ Gzz = 
2j3 that denote the static and time-dependent electric 
field gradients in the respective directions, with the latter 
oscillating at radio- frequency (RF) Vt. In highly excited 
states, ions and atoms cannot be considered point-like in 
such an environment as the trapping field can vary sig- 
nificantly on the size of an electronic orbit as has been 
shown in Refs. [l51[20]. The quantum mechanical treat- 
ment of this problem has been discussed in Ref. [19 but 
in order to be self-contained we briefly summarize the 
main steps: The Hamiltonian of a single ion of mass M 
and charge e reads [10^ 

i7ion(R, r) = i/c(R) + H,(v) + H,,(v) + i7ec(R, r). (1) 

Here i^c(R) = E^[^|/(2M) + MujIrI/2] describes the 
harmonic motion of the center-of-mass (CM) coordinate 
R = (X, y, Z) and the index ^ = X, y, Z refers to 
the respective components. This field results in the 
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FIG. 2. Structural configurations of an electronically excited 
three-ion crystal and a crystal in which all ions are in the elec- 
tronic ground state. The assumed configuration depends on 
the static field gradient /3 and the principal quantum num- 
ber n. While in regions (a) and (c) the configurations are 
identical, an electronic excitation of the central ion (red) is 
accompanied by a structural configuration change in region 

(b). 



trap frequencies ujx = y^if^^ — e(l + £:)/3]/M, ujy = 

^2[^-e{l-e)P]/M and ujz = ^W^. The pa- 
rameter < e < 1 breaks the axial symmetry of the trap 
which simplifies the theoretical discussion as it removes 
some mode degeneracies. Terms describing the coupling 
of the micromotion to the CM are neglected. The Hamil- 
tonian i^e(r) governs the dynamics of the valence electron 
with the relative coordinate r. Its eigenstates 1^^) and 
eigenenergies En (the index n contains all relevant quan- 
tum numbers) are obtained from a model potential calcu- 
lation [21]. Finally, the Hamiltonians i^et(r) = — e<I>(r, t) 
and i^ec(R, r) = — 2eR • G • account for the electron- 
trap interaction and the fundamentally important cou- 
pling between the CM and electronic motion. The to- 
tal Hamiltonian of the three-ion crystal in Fig. [l^ is 
obtained by further including the Coulomb interaction 
among the ions. For the two ground-state ions this is 
straight-forward. In the case of the Rydberg ion the 
Coulomb interaction of the valence electron and the dou- 
bly charged core with the remaining two ions is accounted 
for by a multipole expansions [10 . The dominant contri- 
bution of this expansion is the dipole-charge interaction 
Vdc = [eV(47reo)] E^-(R - ^j) • r/|R - R^f where eo is 
the permittivity of vacuum, R and r denote the CM and 
electronic coordinate of the Rydberg ion and the vectors 
Rj contain the coordinates of the ground state ions. 

When excited to the Rydberg state the coupling 
Hamiltonians Hqc and Vdc create an additional pon- 
deromotive potential for the CM motion of the cen- 
tral ion. This potential provides a tighter radial con- 
finement Vadd = MujI^^{X'^ + Y'^)/2 with the trap 



frequency Wadd = ^/[S{e^a^Y / {Mil'^Y + Ae^a^] x^^f /M. 
Here x^^) = |(M/„|x|vl/ „)| V(i?n - E^) is propor- 

tional to the polarizability of the Rydberg electron. In 
deriving 14dd we made use of the fact that the typical Ke- 
pler orbital frequency of the Rydberg electron is about 
27r X 100 GHz, which is much higher than the RF fre- 
quency Vt (typically 27r x 10... 100 MHz) and also much 
higher than typical trap frequencies (~ 27r x 1 MHz). 
This allows us to calculate the electronic motion within 
a quasistatic approximation where we treat time as a pa- 
rameter [19 . Compared to the vibrational dynamics of 
the ions, however, this surface will oscillate at the fast 
RF frequency. Therefore, the CM of the central ion ex- 
periences an effective ponderomotive potential. 

For a quantitative discussion we choose an RF field 
with gradient a = lO^V/m^ oscillating at a frequency 

= 27r X 30 MHz and a parameter e = 0.15 for splitting 
the degeneracy of the trap frequencies in x and y direc- 
tion. For the following discussion it is more convenient 
to discuss the expected structural configuration change 
in the ion crystal as a function of the gradient P rather 
than in terms of the trap anisotropy A. The critical gra- 
dient /3o at which the linear-zigzag transition occurs for 
a crystal of three ground state ions is given by 
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FIG. 3. (a) Displacement |Ko| of the CM potential surfaces 
corresponding to a ground state and and excited three-ion 
crystal, as illustrated in Fig. [l]:. |Ko| is plotted against the 
trapping field gradient /3 with the abscissa starting at /3o. A 
large displacement affects the Franck-Condon (FC) factors 
that quantify couplings between different vibrational states 
of the two potential surfaces. Panel (b) shows the FC factors 
for transitions from the vibrational ground state of the lower 
potential surface. 
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Due to the additional trapping potential T4dd this criti- 
cal gradient changes when the central ion is excited to a 
Rydberg state. This leads to a state-dependent critical 
gradient Pdn) which is shown Fig. [2] together with Pq. 
The plot can be thought of as a 'phase diagram': in region 
(a) both the excited and the ground state crystal prefer 
a linear arrangement. A zigzag configuration is preferred 
in both cases when gradient and principal quantum num- 
ber are chosen from region (c). In region (b) an electronic 
excitation is accompanied by a structural configuration 
change. 

In the remaining part of this work, we will analyze the 
spectroscopic properties of the crystal which derive from 
this state dependent configuration change. We proceed 
by introducing normal modes that describe the vibra- 
tional motion of the ion crystal on the two potential sur- 
faces which correspond to the electronically excited ion 
crystal and the crystal with all ions in the ground state 
(see sketch in Fig. [l]^). This results in the vibrational 
Hamiltonian i^vib = ^^fux)ia\aiP- + ^j^^j^l^jP- 
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Here ai{bj) are the normal mode annihilation operators 
on the lower (upper) potential surface with vibrational 



frequencies uJi {uj 



(1 ^ crz)/2 projects on the 



electronic ground state/Rydberg state of the central ion 
and involves the Pauli matrix . To simplify the discus- 
sion we assume that the central ion is already excited to 
the metastable 31^3/2 state, which can be done without 
affecting the crystal configuration. After including the 
Rydberg excitation laser on the 3I)3/2-'^^i/2-transition 



with Rabi frequency Qj, the total Hamiltonian becomes 
H = i^vib + where 



h.c. 



(2) 



and a± = l/2{ax ± i(Jy) (with 1 31)3/2) = \i) and 
1^^1/2) = It))- Here we have introduced the detuning 
hS = huji, — {Enp^^^ + AE + Eqi) where cjl is the laser 
frequency, EnP^^^ is the nPi/2 state energy relative to 
the metastable 3D3/2 state, AE is the energy shift of 
the two potential surfaces resulting from the structural 
change of the ions, and ^et a constant energy shift due 
to the electron-trap coupling i^et- The coefficients 
are the so-called Frank- Condon (FC) factors and given 
by the overlap integrals of the vibrational modes of the 
lower and excited potential surface. Here [n] and [m] are 
multi-indices containing all phonon occupation numbers 
of the lower/upper modes. They are also used for ab- 
breviating a^n] = Yin ^nii whcrc is the annihilation 
operator corresponding to the n-phonon state of the z-th 
mode. The same notation applies to 

The FC factors crucially depend on the mutual dis- 
placement |Ko| of the minima of potential surfaces. In 
the region (b) of the 'phase diagram' in Fig. |2] where the 
configurations of the ground state and the excited crys- 
tal differ substantially, |Ko| can significantly exceed the 
harmonic oscillator length, typically ~ 10 nm. This is 
illustrated in Fig. [Sk for the Rydberg 38Pi/2 state. In 
Fig. we show the corresponding FC factors for tran- 
sitions starting from the vibrational ground state of the 
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FIG. 4. Due to non-zero FC factors phonon states contain- 
ing many vibrational quanta can be directly addressed. The 
data in the main figure shows the excitation probability of 
an 11-phonon state in the upper potential surface. We have 
considered the 38Pi/2 state of ^°Ca^, the laser is resonant 
with the desired phonon state and Ql = 27r x 0.2 MHz, 
13 = 1.1392 X lO'^V/m^ The corresponding FC factor is 
^ 0.195. The inset (a) shows the phonon energies of the low- 
est energy normal modes, cjq (all ions in the electronic ground 
state) and (electronically excited crystal). The inset (b) 
shows the density profile (cut along ^ = 0) of the Rydberg 
ion CM motion in the excited vibrational 11-phonon state. 



lower potential surface, which are numerically calculated 
by a Gaussian integral method [22). We observe FC fac- 
tors that indicate significant overlap with highly excited 
vibrational states of the upper potential surface. The 
calculation furthermore shows that only the lowest vibra- 
tional mode in each potential surface determines the FC 
factor. Higher modes are virtually identical and the FC 
factors for transitions changing their vibrational state are 
zero. The frequencies of the two lowest modes {ujq and 
uo) are shown in the inset (a) of Fig. [4j 

The FC factors in Fig. ^jp are a result of strong me- 
chanical forces that accompany the Rydberg excitation 
of the ion crystal. These forces are substantial also for 
crystals with large vibrational frequencies ~ MHz. This 
is a striking difference to the common situation encoun- 
tered in ion traps where merely low-lying electronic states 
are laser excited. Here one is usually in the Lamb-Dicke 
regime and the direct excitation of high-lying vibrational 
states is strongly suppressed. We can use the emerging 
large mechanical forces to selectively populate high-lying 
vibrational state of the ion crystal. This is depicted in 
Fig. [4] for a state on the upper potential surface con- 
taining 11 vibrational quanta which is addressed directly 
from the vibrational ground state of the lower surface. 
The process is state-selective as long as the product of 
the FC factor and the excitation Rabi frequency Qj, is 
smaller than the phonon energy ujq. Here we chose a re- 
alistic Rabi frequency of (^l = 27r x 0.2 MHz [19 which 



achieves a population of the 11-photon state with more 
than 97% probability in t ^ 13 /j^s. This process creates 
a 'moving' Rydberg ion with a probability density shown 
in Fig. [4]3. 

This situation is distinct from conventional cold neu- 
tral Rydberg systems which are often studied in the so- 
called 'frozen' regime [23l [24] and where motion even- 
tually leads to decoherence caused by Penning ioniza- 
tion [25 . In our case the 'vibronic dynamics' is coherent 
due to the absence of collisions and the long lifetime of 
^^Ca+ ions in the nPi/2-state. For n > 15 this lifetime 
is To = 1.534 X (n — (^i)^"^^^ns, with the quantum defect 
Si ^ 1.436. These numbers agree well with qualitative 
estimates found in Ref. [2F and yield tq ~ 66 /is for a 
38Pi/2 state. 

In conclusion, the excitation of a cold ion crystal to 
high-lying electronic states was studied and shown to 
be accompanied by large coherent forces. This can find 
practical application in the realization of geometric quan- 
tum gates [27 and the coherent separation of long ion 
chains. Excited cold ion crystals permit moreover the 
quantum simulation of coherent molecular phenomena 
and in principle also allow the inclusion of radiative tran- 
sitions between 'molecular' energy surfaces as tight traps 
can grant confinement even after the de-excitation of Ry- 
dberg states by incoherent photon emission. 
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